Anisotropic Monopole Heat Transport in the Spin-Ice Compound Dy2Ti207 
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We report a study of the thermal conductivity k of the spin-ice material Dy2Ti207. From the 
anisotropic magnetic-field dependence of k, and by additional measurements on the phononic ref- 
erence compounds Y2Ti207 and (Dyo.5Yo.5)2Ti207, we are able to separate the phononic and the 
magnetic contributions to the total heat transport, i.e. Kph and Kmag, respectively, which both 
depend on the magnetic field. The field dependent Kph(-B) arises from lattice distortions due to 
magnetic-field induced torques on the non-coUinear magnetic moments of the Dy ions. For Kmag, 
we observe a highly anisotropic magnetic-field dependence, which correlates with the correspond- 
ing magnetization data reflecting the different magnetic-field induced spin-ice ground states. The 
magnitude of ^mag increases with the degree of the ground-state degeneracy. This anisotropic field 
dependence as well as various hysteresis effects suggest that Kmag is essentially determined by the 
mobility of the magnetic monopole excitations in spin ice. 

PACS numbers: 66.70.-f, 75.40.Gb 75.47.-m 



I. INTRODUCTION 

Dy2Ti207 is a geometrically frustrated spin system 
with a degenerate, so-called spin-ice ground state. The 
magnetic Dy sites in Dy2Ti2 07 form a pyrochlore lattice, 
consisting of corner-sharing tetrahedra. A strong crys- 
tal field results in an Ising anisotropy of the magnetic 
moments of the Dy ions, which align along their local 
easy axes in one of the {111} directions and point either 
into or out of the tetrahedra. Possible ground states in 
zero magnetic field are given by the "ice rule" : two spins 
point into and two out of a tetrahedron. This behavior is 
analogous to the hydrogen displacement in water ice re- 
vealing a residual zero-temperature entropj^^^— . Magnetic 
excitations can be created pairwise by flipping one spin 
resulting in two neighboring tetrahedra with configura- 
tions lin-3out and 3in-lout, respectively. Such a dipole 
excitation can fractionalize into two individual monopole 
excitations which can freely propagate in zero magnetic 
field^-ii. 

Recently, we showed that these exotic magnetic exci- 
tations contribute to the heat transport in Dy2Ti207^. 
At temperatures around ~ 0.6 K, the magnetic contri- 
bution Kmag accounts for almost 50% of the total ther- 
mal conductivity k. A magnetic field i3||[001] lifts 
the ground-state degeneracy of the magnetic system in 
Dy2Ti207. This results in a complete suppression of 
Kmag for rather small magnetic fields of ~ 0.5 T. An 
alternative interpretation had been proposed earlier in 
Ref. [TqI . where the thermal conductivity was assumed to 
be purely phononic and the field dependence of k was 
attributed to phonon scattering on magnetic excitations. 
Our data, however, do not support this interpretation 
(c/. discussion in Ref. [H). In Ref. [H, we studied Kmag 
for a magnetic field applied parallel to [001]. By apply- 
ing the magnetic field parallel to [111] or [110], one can 
realize more complex field-induced ground states with 
different degrees of degeneracy. For B\ \ [111], the field- 
induced ground state below 1 T (Kagome-ice state) is 



3-fold degenerate consisting of different types of 2in-2out 
tetrahedra— Above 1 T, the ground state changes 
into a non-degenerate ground state consisting of alter- 
nating 3in-lout and lin-3out tetrahedra. For a magnetic 
field parallel to [110], two spins per tetrahedron are per- 
pendicular to B and are, thus, not affected. This leads 
to a 2-fold degenerate ground state, even for arbitrarily 
large magnetic fields. In this paper, we study the de- 
pendence of the magnetic heat transport on the degree 
of the field-induced ground-state degeneracy realized by 
the field directions [001], [111], and [110]. 

In the course of this study, we also extend the 
field-dependent thermal-conductivity measurements of 
Dy2Ti207 to magnetic fields up to 7 T. An addi- 
tional field dependence of n is observed for higher 
fields (depending on the actual field direction) which 
we identify as a field-dependent phononic background 
Kph(-B) of Dy2Ti207. The phononic background is essen- 
tially given by the field-dependent k{B) of the Y-doped 
(Dyo.5Yo.5)2Ti207. This compound can be regarded as a 
magnetic reference system with strongly suppressed spin- 
ice features. Compared to the non-magnetic Y2Ti207, 
the half-doped reference compound can be utilized to 
study the magnetic-field dependence of the phononic con- 
tribution Kph of Dy2Ti207. The field dependence of Kph 
can be explained by a magnetostrictive lattice distortion 
originating from torques of the Dy momenta induced by 
the external magnetic field. 



II. EXPERIMENTAL 

Single crystals were grown from sintered bars of Ti02 
(3N, Sigma-Aldrich) and Y2O3 (4N, Alfa Aesar) and/or 
Dy203 (4N, REacton) in proper stoichiometry. The 
sintering processes were done in air at 1400 °C over 
night using corundum boats as crucibles. For the crys- 
tal growth, a fioating-zone technique was applied inside 
a four-mirror image furnace to acquire centimeter-sized 
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single crystals. To obtain crystals without cracks, the 
floating zone was run twice trough the bars: first at a 
rate of 25mm/h and subsequently with 7nim/h. Ambi- 
ent pressure of pure oxygen was chosen as atmosphere 
during the growth. The pale yellow crystals exhibit 
high reflectance at a grazing angle but otherwise full 
transparency. Dy2Ti207 is slightly more yellow than 
(Dyo.5Yo.5)2Ti207 and YzTiaOv. The resulting crys- 
tals were found to be single domain by extensive Laue- 
photo investigations and the bulk purity was checked by 
X-ray powder diffraction using a STOE D5000 diffrac- 
tometer, Cu Kq, radiation, and Bragg-Brentano reflec- 
tion mode. The correct composition was confirmed by 
energy-dispersive spectroscopy inside a scanning electron 
microscope. 

The thermal conductivity was measured using the 
standard steady-state method. The temperature differ- 
ence was produced by a heater attached at one end of the 
sample and measured by a pair of matched Ru02 ther- 
mometers. The thermal-conductivity measurements were 
performed on bar-shaped single crystals of approximately 
3x1x1 mm'^ with the long edge parallel to the [110] di- 
rection. The heat current was directed along the long 
edge of the crystal. In most cases, the magnetic fields 
applied along the different directions were oriented per- 
pendicular to this long edge. Hence, considerable de- 
magnetization effects have to be taken into account. The 
demagnetization field is calculated on the basis of exper- 
imental magnetization data which were measured with 
a home-built Faraday magnetometer on thin samples to 
minimize demagnetization effects within the magnetiza- 
tion measurements. The magnetostriction measurements 
were done with a home-built capacitance dilatometer, 
which allows to measure the uniaxial length changes ei- 
ther parallel or perpendicular to the applied magnetic 
field. 



III. RESULTS 

Fig.[l]shows the relative change of the thermal conduc- 
tivity k{B)/k{OT) of Dy2Ti207 for B\\ [111] at temper- 
atures between 0.35 K and 1 K. The different datasets 
are shifted with respect to each other. With increasing 
field, k{B) decreases and exhibits a plateau within the 
Kagome-ice phase below 1 T. For higher fields, k{B) 
further decreases and shows a kink around 1.5 T. The 
plateau feature is more pronounced at lowest tempera- 
ture and broadens towards higher temperature. The inset 
of Fig. [T] displays the magnetization M{B) of Dy2Ti207 
for B\ \ [111] at 0.4 K, which is in agreement with pre- 
vious studies, e.g. Refs. 4 and [2l|. The magnetization 
as well as the thermal conductivity have been measured 
with slow sweep rates^^ to avoid thermal-runaway effects, 
c/. Ref. [Tl|. The basic features of the k{B) data are re- 
fiected within the magnetization data, which also exhibit 
a pronounced plateau within the Kagome-ice phase. At 
1 T, M{B) shows a sharp kink and increases towards the 



T ' r 




H„(H-DM) (T) 

FIG. 1. (Color online) Relative change k/ko of Dy2Ti207 
for _B II [111] at various constant temperatures. The different 
curves are shifted by 0.15 with respect to each other. Inset: 
Magnetization M{B) of Dy2Ti207 at 0.4 K. The arrows in- 
dicate the field-sweep direction. Demagnetization effects are 
taken into account. 



saturation value, which is essentially reached at ^ 1.5 T. 
This increase of M{B) is accompanied by a decrease of 
K between 1 T and 1.5 T. The hysteresis in M(B) at 
low temperature (~0.4K) below '^0.3T can also be 
observed in the thermal-conductivity data. Even the 
remnant zero-field magnetization observed in the mea- 
surement with decreasing field (open symbols), which 
results from slow relaxation processes of the spin ice, 
is refiected within the thermal-conductivity data. Be- 
low ~ 0.4 K, the K data for decreasing field end in a re- 
duced zero-field value, which slowly relaxes back to the 
initial zero-field value (after zero-field cooling). Such a 
correlation of magnetization and thermal conductivity 
has already been observed for B \ \ [001] (Ref. jl^). For 
B II [111], however, clear differences between M{B) and 
k{B) are also found. First of all, at lowest temperature, 
k{B) within the Kagome-ice phase strongly depends on 
whether the measurement is performed with increasing 
field (after cooling in zero field) or with decreasing field 
(starting from high magnetic fields) . Such a hysteresis is 
not present in the magnetization data. Even at lowest 
temperature, the magnetization within the Kagome-ice 
phase shows no hysteresis with respect to the field-sweep 
direction. Secondly, M{B) and n{B) show different high- 
field behaviors. At 0.4 K, the magnetization is almost 
saturated above 1.5 T, whereas the thermal conductivity 
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H„(H-DM) (T) 

FIG. 2. (Color online) Relative change k/kq of Dy2Ti207 
for B II [001] at various constant temperatures between 0.4 K 
and 8 K. Demagnetization effects are taken into account. 
Panel (a) compares the field dependencies of magnetization 
and thermal conductivity at 0.4 K. The complete suppression 
of Kmag below 1 T is tagged by the dashed lines. The inset is 
an enlargement of the low-field region, where the field-sweep 
directions are marked by arrows. 



k(B) further decreases with increasing field even above 
1.5 T. The question arises whether these effects can be 
attributed to the phononic or to the magnetic contribu- 
tion of K. 

In order to analyze this, we will first discuss the high- 
field dependence of k for magnetic fields above ~ 1.5 T. 
Because the magnetization is essentially saturated above 
1.5 T, the magnetic system in Dy2Ti207 is fully polar- 
ized. Hence, one would expect the field dependence of 
K above 1.5 T to originate from a change of Kph rather 
than from Kmag- By means of the data shown in Fig. [TJ 
however, one cannot straightforwardly separate Kmag and 
Kph. In Ref.[l^ we showed that Kmag vanishes for a rather 
small field of ~ 0.5 T applied parallel to [001]. Compared 
to i? II [111], the field-induced ground state for B\ \ [001] 
is less complex, as the ground-state degeneracy is lifted 
in one single step. Therefore, we start with the anal- 
ysis of the high-field thermal conductivity of Dy2Ti2 07 
for B II [001]. The data published in Ref. [H were only 
measured up to a (demagnetization-corrected) maximum 
field of 0.5 T. Here, we extend the present k data to mag- 
netic fields up to 7 T. The results are displayed in Fig. [21 
where the relative change k{B)/k{OT) is shown for tem- 
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FIG. 3. (Color online) Comparison of the zero-field thermal 
conductivity of Dy2Ti2 07 (closed blue circles) with the non- 
magnetic reference Y2Ti207 (open red squares) and the half- 
doped (Dyo,5Yo,5)2Ti207 (open green triangles). 



peratures between 0.4 K and 8 K. Below 0.5 T, the data 
shown in Fig. [5] confirm the data previously published 
in Ref. [H. Fig. [Sfa) directly compares the field-induced 
change of k with the magnetization M{B) at the lowest 
measured temperature of 0.4 K. Below 0.5 T, k{B) step- 
likely decreases. As can be seen in the inset of panel (a), 
the sharp drop of k{B) below 0.5 T clearly correlates 
with M{B) and can, thus, be attributed to Kmag- Above 
^ 1.5 T, k{B) further decreases. This additional field de- 
pendence k{B), however, cannot be attributed to Kmag as 
all spins are fully polarized. This is shown exemplarily 
for 0.4 K, where the magnetization is essentially satu- 
rated above 1.5 T (apart from a slight linear increase, 
see below). The field dependence of k also cannot be 
explained by phonon scattering on magnetic excitations 
as the excitation gap monotonically increases with an in- 
creasing field B II [001]. Fig.^b) shows the field-induced 
relative change of k{B) at higher temperatures up to 8 K. 
Above ~ 4K, Kmag vanishes, whereas a distinct high-field 
dependence of k above 1 T is still observed even at 8 K. 
As will be shown below, this high-field dependence of the 
thermal conductivity of Dy2Ti207 can be attributed to a 
magnetic-field dependent phononic contribution Kp\^{B). 

In our previous worki^, we studied the non-magnetic 
Y2Ti207 as a purely phononic reference compound to 
the magnetic spin ice Dy2Ti207. At temperatures below 
^ 3 K, the thermal conductivity of ¥2X1207 shows a very 
similar power-law behavior as the phononic background 
of Dy2Ti207 which was obtained by applying a mag- 
netic field of 0.5 T parallel to [001]. For higher tempera- 
tures, however, the thermal conductivity of Dy2Ti207 
is significantly smaller (see Fig. [3]). Most likely, this 
suppression originates from additional phonon scatter- 
ing on the crystal-field excitations of the 4f electrons of 
the Dy ions. Here, we present an alternative approach 
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in order to study the phonon background of k by using 
the half-doped compound (Dyo.5Yo.5)2Ti207. As every 
second Dy ion is replaced by a non-magnetic Y ion, i.e. 
on average two of the four Dy ions per tetrahedron are 
replaced by Y ions, the spin-ice properties are supposed 
to be strongly suppressed. However, due to the similar 
ionic radii of Dy'^"'' and Y^"^ , one can assume both com- 
pounds to have comparable phononic properties. Fig. [3] 
shows the zero-field thermal conductivity of the spin-ice 
compound Dy2Ti207 together with the thermal conduc- 
tivity of the non-magnetic reference Y2Ti207 and the 
half-doped (Dyo.5Yo.5)2Ti207. In all measurements, the 
heat current has been directed along the [110] direc- 
tion. Below 1 K, both phononic reference compounds 
have very similar k values well below k of Dy2Ti2 07. At 
higher temperature, the half-doped compound shows a 
behavior similar to the spin-ice compound, rather than to 
the non-magnetic Y2Ti207, which has significantly larger 
K values. This supports our interpretation of an addi- 
tional phonon scattering on the crystal-field excitations 
of the Dy ions, as given above. This additional scat- 
tering, indeed, can also be observed for the half-doped 
(Dyo.5Yo.5)2Ti207, leading to k values close to the val- 
ues of Dy2Ti207. The slightly smaller values of the half- 
doped compound can be explained by an enhanced defect 
scattering of phonons which originates from the partial 
Dy substitution. 



IV. DISCUSSION 
A. Phononic Heat Conductivity Kph 

The main advantage of the half-doped compound 
(Dyo.5Yo.5)2Ti207, compared to Y2Ti207, is the fact 
that it is still strongly magnetic and can, therefore, be 
utilized to study the magnetic-field dependence of the 
phononic background Kph(-B) of Dy2Ti207. As shown 
exemplarily for B\ \ [111] in Fig.lHd), k of Y2Ti207 (red 
squares) has no magnetic-field dependence, as expected 
for a non-magnetic compound. Figs. l3Ua)-(d) compare 
the thermal conductivity k{B) of Dy2Ti207 (blue circles) 
and of (Dyo.5Yo.5)2Ti207 (green triangles) for various 
magnetic-field directions, where the heat current j is al- 
ways directed parallel to [110]. First of all, even at lowest 
temperature, k{B) of (Dyo.5Yo.5)2Ti207 is not hysteretic 
with respect to the field-sweep direction (not shown). 
This is contrary to the undoped Dy2Ti207, where at low- 
est temperature 0.4 K), a clear hysteresis in k{B) is 
observed below ~ 0.3 T for all considered field directions, 
see e.g. Figs. [1] and [5j 

The data in Figs. Hl^a) and (b) are measured with 
a magnetic field parallel to [001]. At 2 K [Fig. gja)], 
k{B) of Dy2Ti207 shows two steplike anomalies at 
~0.5T and ~2T and continuously decreases above 
the second anomaly. The step at ~ 2 T and the 
subsequent continuous decrease are very well refiected 
by the k{B) data of (Dyo.5Yo.5)2Ti207. The first 
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FIG. 4. (Color online) (a)-(d): Field-dependent k{B) of 
Dy2Ti207 (blue circles) and of the half-doped phononic refer- 
ence compound (Dyo.5Yo.5)2Ti207 (green triangles) for differ- 
ent magnetic-field directions; panel (d) also contains k(_B) of 
Y2Ti207. The inset in panel (b) compares «;(-B) with M{B) 
for (Dyo.5Yo.5)2Ti207. (e) and (f): Field-induced length 
changes of Dy2Ti207 for B [j [001] and B || [111], respectively. 
The inset in panel (e) shows the corresponding field derivative 
of AL/L, where the anomaly at 2 T (arrows) is better seen. 
In all cases, demagnetization effects are taken into account 
and, for clarity, only data obtained with increasing magnetic 
field are shown. 



anomaly (at ^0.5T), which we attribute to the sup- 
pression of Kniag in Dy2Ti207, is, however, absent in 
the (Dyo.5Yo.5)2Ti207 data. Fig. [3Kb) shows the data 
for the same field direction B \ \ [001] at a lower tempera- 
ture of 0.4 K. Here again, the high-field behavior is very 
similar for both compounds, whereas the sharp decrease 
below ^ 0.5 T is absent in the (Dyo.5Yo.5)2Ti207 data. 
The inset of Fig. [U^b) illustrates the different origins of 
the field dependencies of the half-doped compared to the 
undoped compound. One can very clearly see that for 
(Dyo.5Yo.5)2Ti207, the change of k.{B) does not corre- 
late with the magnetization M{B), which is essentially 
saturated at ~ 1 T, whereas k{B) shows a distinct field 
dependence up to ^ 3 T. This is different from the case 
of Dy2Ti207, where the change of M (B) at low field (be- 
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low ^ 0.5 T) directly correlates with the suppression of 
^^mag [inset of Fig. ^a)]. Thus, we conclude that the 
high-field decrease of k of Dy2Ti207 can be identified as 
a field-dependent phononic background Kph{B). 

It is notable, however, that the zero-field k values of 
(Dyo.5Yo.5)2Ti207 at lowest temperature are larger than 
the phononic background of Dy2Ti207 obtained from the 
plateau value around 1 T [Fig. Hljb)]. Most likely, this 
originates from a small remnant magnetic contribution 
^mag in (Dyo.5Yo.5)2Ti207, which is more pronounced at 
lower temperature (0.4 K) than at higher temperature 
[2 K, Fig. Ilja)]. Hence, one cannot directly identify the 
low-field (below ^ 0.5 T) k values of the half-doped com- 
pound as the phononic background Kph of the undoped 
mother compound. Nevertheless, the remnant magnetic 
contribution in (Dyo.5Yo.5)2Ti207 is strongly suppressed 
compared to Kmag in Dy2Ti207, as expected for half dop- 
ing. 

In Fig. IHc), the field dependence of k is shown for 
a magnetic field parallel to [110]. For this particular 
field direction, two of the four spins per tetrahedron in 
Dy2Ti207 are perpendicular to B and are, thus, not af- 
fected. This results in a 2-fold degenerate ground state 
for this field direction and we find that the decline of k{B) 
below 0.5 T is less pronounced as compared to i? || [GDI] 
[dashed line, cf. panel (b)]. This indicates that Kmag can- 
not be completely suppressed by a field applied parallel 
to [110]. Furthermore, above ~0.5T, k is hardly field 
dependent. This is different from the other considered 
field directions, where k{B) shows a distinct high-field 
dependence of k. However, here also, the thermal con- 
ductivity of the half-doped reference compound closely 
resembles the high-field data of Dy2Ti207 at least up 
to 4 T, whereas there are some deviations appearing at 
higher field. The data in Fig. SJ^c) might suggest that 
the zero-field k value of (Dyo.5Yo.5)2Ti207 essentially re- 
fiects the phononic background for Dy2Ti2 07. This is, 
however, misleading as both compounds exhibit a non- 
vanishing zero-field magnetic contribution to k (as shown 
above for B \ \ [GDI]) and, obviously, the zero-field k is in- 
dependent on the magnetic-field direction. Consequently, 
the field-induced 2-fold degenerate ground state realized 
by i? II [IIG] leads to a considerable magnetic contribu- 
tion Kmag, which is hardly field dependent up to highest 
fields of 7 T. 

More complex ground states are realized by a mag- 
netic field applied parallel to [111]. Within the Kagome- 
ice phase below 1 T, the ground state is 3-fold degener- 
ate. Above the transition at 1 T, the magnetic system 
changes to a non-degenerate, fully polarized ground state. 
The highly degenerate Kagome-ice state results in a pro- 
nounced plateau within the field-dependent change of k 
observed for Dy2Ti207 at lowest temperature (Fig. [T]). 
The high-field dependence of k for i3||[lll] is shown 
in Fig. mjd) together with the field-dependent data for 
the half-doped (Dyo.5Yo.5)2Ti207 and the non-magnetic 
Y2Ti207. As mentioned above, Y2Ti207 can only be 



used to estimate the zero-field phononic contribution, as 
K of Y2Ti207 shows no field dependence. In complete 
analogy to the field directions discussed above, the field 
dependence of k above 2 T is very well described by n of 
the half-doped reference compound, whereas the spin-ice 
features at lower field are almost completely suppressed. 
This, again, confirms our interpretation of k(B) of the 
half-doped reference compound as the field-dependent 
phononic background of Dy2Ti207. However, it is more 
difficult to distinguish between the magnetic and the 
phononic contributions of Dy2Ti207 for B\ \ [111]. The 
field dependence well above ^ 1.5T can certainly be at- 
tributed to Kph, but around 1.5 T, Kph(i?) and the sup- 
pression of Kmag due to the lifting of the ground-state 
degeneracy above the Kagome-ice phase overlap. Con- 
sequently, the suppression of Kmag appears as a kink in 
k{B) for i3|][lll], rather than a pronounced plateau, 
as observed for B\ \ [001]. The kink is located close to 
the zero-field phononic background obtained from the 
B 11 [001] data [horizontal dashed line, cf. Fig. i^b)]. 
This suggests that the field dependence Kph{B) below 
^ 1.5T is much smaller than the suppression of Kmag- 
Hence, when we restrict the discussion to the low-field 
region below 1.5 T, we can assume a constant, i.e. field- 
independent, phononic background Kph of Dy2Ti207. 
The magnetic contributions Kmag extracted on the basis 
of this assumption are discussed in the following subsec- 
tion. 

A possible explanation of the field dependence of the 
phononic background Kph of Dy2Ti207 is given in terms 
of field-induced lattice distortions which originate from 
torques jl x B of the Dy momenta jl induced by an ex- 
ternal magnetic field B. As the different local easy axes 
are not collinear, these torques emerge for any field di- 
rection. The resulting lattice distortions can be visual- 
ized by magnetostriction measurements. Figs. SKe) and 
(f) show the field-induced length changes of Dy2Ti207 
for B\ \ [001] and B \ \ [111], respectively. In both cases, 
the Dy2Ti207 crystals elongate parallel to the magnetic- 
field direction (red squares), where the effect is more pro- 
nounced for B\ \ [001]. In the high-field region, i.e. for 
magnetic fields above the (almost) saturation of M{B), 
this elongation is essentially linear in B and temperature 
independent in the temperature range between G.25 and 
2 K (not shown). Thus, the high-field magnetostriction 
cannot be attributed to the spin-ice system in Dy2Ti207. 
For i? II [111] [Fig. IHf)], we also measured the field- 
induced length changes perpendicular to B, namely par- 
allel to [IIG] and [112], and we find that the crystal con- 
tracts in both perpendicular directions, but to a much 
smaller extent as compared to the elongation parallel to 
B\ \ [111]. Thus, the magnetic field causes a strong distor- 
tion of the lattice due to torques of the large Dy momenta 
which tend to align parallel to the external field. These 
lattice distortions influence the phononic properties, in 
particular, the phononic thermal conductivity. The step- 
like anomaly observed in k{B) for B\ \ [001] around 2 T 
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[Figs, mja) and (b)] has, indeed, an equivalent within the 
magnetostriction data, which show a smah anomaly at 
the same magnetic field. In Fig. |4lje), an arrow marks 
this anomaly, which is better seen in the field-derivative 
of the length change [inset of Fig. IH^e)] . The origin of this 
anomaly is, however, unclear. In sum, the field depen- 
dence of the thermal conductivity of Dy2Ti2 07 observed 
in the high-field region [Figs, ma)-(d)] can be explained 
by lattice distortions and can, thus, be attributed to a 
field-dependent phononic contribution Kph(i?). In addi- 
tion, the weak, essentially linear high-field increase of the 
magnetization observed above about 1.5 T [c/. Fig.[5][a)] 
can also be explained by the torques fix B. For example, 
the additional increase of M{B) between 2 and 7 T can 
be traced back to a tilting of the Dy momenta by ^ 1.5°. 



B. Magnetic Heat Conductivity 



Now we proceed with the discussion of the magnetic 
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contribution 



^mag 



of Dy2Ti207. In the previous para- 



graphs, we showed that 



can be suppressed by a 



rather small magnetic field, depending on the actual field 
direction. As Kph is only weakly field dependent for 
small magnetic fields, we can extract Kmag by assuming 
a constant, i.e. field- independent, phononic background 
[dashed line, c/. Fig. HJb)]. The field-dependent n 



mag 



for 

the different field directions are displayed in Figs, 
and (b) at 1 K and 0.4 K, respectively. For clarity, we 
only show the data measured with decreasing magnetic 
field, which are (for B > Q) closer to thermal equilib- 
rium than the data obtained with increasing field (see 
discussion below). The magnetic contribution is maxi- 
mum in zero field (6-fold degenerate ground state) and 
completely vanishes for B\ \ [001] (non-degenerate). For 
B II [110] (2-fold degenerate), Kmag cannot be completely 
suppressed, even for a rather large field of ~ 7T. For 
this particular field direction, the crystallographic direc- 
tions [110] and [110] become inequivalent. The spins 
with a non-vanishing component parallel to B become 
fully polarized and form so-called a chains running par- 
allel to B. The remaining spins are perpendicular to B 
and form /3 chains that are oriented perpendicular to B, 
i.e. a\\ [110] and /3|| [110], c/. Refs.llltll and [13. For 
-B|| [110], we measured Kmag within two different mea- 
surement setups, either with the heat-current direction 
along the a chains, i.e. j\\B^ or with the heat current 
along the (5 chains, i.e. j _L _B. At 1 K, both datasets 
agree with each other [Fig. [51[a)], whereas at 0.4 K, K,„ag 
is significantly larger along the /? chains [Fig. [Sjb)]. A 
direct comparison of the field-induced relative changes 
of K for j 1 1 a and j 1 1 /3 at IK and 0.4 K is displayed 
in Figs, ^c) and (d), respectively. We note that the 
data of Figs. [5jc) and (d) have been obtained in a dif- 
ferent measurement run on another crystal as compared 
to the data displayed in panels (a) and (b). The low- 
temperature (0.4 K) data for both heat-current direc- 
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FIG. 5. (Color online) (a) and (b): Field-dependence of Kmag 
of Dy2Ti207 for different magnetic-field directions at 1 K and 
0.4 K, respectively. For clarity, only the data measured with 
decreasing field are shown, (c) and (d): Field- induced relative 
change fi:(_B)/ft:(0 T) for i3|| [110] and different directions of 
the heat current j at 1 K and 0.4 K, respectively. The arrows 
indicate the different field-sweep directions. Demagnetization 
effects are taken into account. 



tions show a pronounced hysteresis with respect to the 
field-sweep direction and the k{B — 0) data measured 
with decreasing field end in a reduced zero-field k value, 
as it is also observed for B\ \ [001] and B \ \ [111]. More- 
over, the low-temperature data for B \ \ [110] clearly reveal 



that 



''mag 



is larger for a heat current along the /? chains 



than along the a chains. Such an anisotropy of Kmag 
with respect to j can be qualitatively explained within 
the microscopic model of monopole excitations propa- 
gating via single spin flips. The spins of the a chains 
are easily polarized by i3|| [110], whereas the spins of 
the /3 chains are perpendicular to B and are, thus, not 
affected. As a consequence, the monopole propagation 
along the a chains should be suppressed compared to the 
/3 chains. Within the model of single spin flips, one would 
even expect that monopole excitations can only propa- 
gate along the fj chains, which are separated by (fully 
polarized) a chains, whereas Kmag should completely van- 
ish for j II a. A complete suppression of Kmag along the 
a chains is, however, not observed [Figs. [5l[a) and (b)]. 
This shows that the single spin-flip formalism, which only 
accounts for nearest-neighbor interaction, is an oversim- 
plification and the long-range dipolar interaction has to 
be taken into account. 

For a magnetic field parallel to [111] within the 
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Kagome-ice phase below 1 T (3-fold degenerate), Kmag 
is significantly larger than for B\ \ [110] [Figs, ^a.) and 
(b)]. Above 1 T, the ground-state degeneracy is lifted, 
resulting in a suppression of Kmag, which essentially van- 
ishes above 1.5T (indicated by a kink). As men- 
tioned above, the field dependence of Kph cannot be 
neglected above ~ 1.5 T. Hence, the assumption of a 
field-independent phononic background is no longer jus- 
tified, and Kniag extracted here does not exactly ap- 
proach zero. The hysteretic behavior below 1 T, how- 
ever, can certainly be attributed to Kmag because no 
hysteresis of Kph(-B) is observed for the reference com- 
pound (Dyo.5Yo.5)2Ti207. At lowest temperature (below 
^ 0.4 K), Kmag within the Kagome-ice phase (below 1 T) 
depends on whether the measurement is performed with 
increasing field after zero-field cooling or with decreasing 
field starting from the fully polarized phase (above 1 T) 
(c/. Fig.[T]). The hysteresis in k{B) below ~ 0.3 T and the 
fact that the measurement with decreasing field does not 
end at the initial zero-field-cooled value is also present 
for B II [110] [Fig. Eld)] and for i?|| [001] (c/. Ref. HI), 
and there are corresponding hysteresis effects in the re- 
spective magnetization data (inset of Fig. [T]). 

The hysteresis within the Kagome-ice phase, however, 
has no equivalent within the magnetization data. The 
absence of such a hysteresis in M{B) can be explained 
by the single-tetrahedron model. The 3-fold degener- 
ate ground state consists of three different, energetically 
equivalent tetrahedron configurations. As these configu- 
rations have the same magnetization component parallel 
to B II [111], they are indistinguishable for this partic- 
ular field direction. The thermal conductivity, however, 
seems to be sensitive to the actually realized ground-state 
configuration. Above 1 T, the spin ice is fully polarized, 
and the thermal conductivity does not depend on the ac- 
tual field-sweep direction. When entering the Kagome- 
ice phase from higher field, the spins of the triangular 
planes which separate the Kagome planes stay polarized, 
whereas the spins of the Kagome planes arrange in such 
a way that the tetrahedra obey the ice rule 2in-2out. As 
this configuration can easily be achieved by flipping only 
one spin per tetrahedron, it can be assumed that this 
configuration is close to thermal equilibrium. As a con- 
sequence, Kmag stays almost constant on the Kagome-ice 
plateau. The situation is more complicated when ap- 
proaching the Kagome-ice phase from below, i.e. from 
the highly entropic ground state after zero-field cooling. 
The magnetic field i?|| [111] initially polarizes the spins 
of the triangular planes. Flipping such a spin, however, 
produces a monopole/anti-monopole pair in two neigh- 
boring Kagome planes. To reach thermal equilibrium, 
these monopoles (and anti-monopoles) have to be annihi- 
lated. This out-of-equilibrium state results in a reduced 
Kmag) and time-dependent measurements of Kmag (not 
shown) reveal that it very slowly converges towards the 
larger equilibrium value. When approaching the transi- 
tion at 1 T, the energy gap to the fully polarized state 
vanishes and, as a consequence, the Kagome planes can 



easily reach thermal equilibrium. This results in an in- 
crease of K,{B) towards 1 T. An interesting finding is 
the fact that in the field- decreasing run, i.e. in (or at 
least close to) thermal equilibrium, Kmag remains con- 
stant on the Kagome-ice plateau although the energy 
gap to the (anti-)monopole excitations increases with de- 
creasing field. 

In sum, our data for the different magnetic-field direc- 
tions reveal that there is a magnetic contribution Kmag, 
whose magnitude directly correlates with the degree of 
the ground-state degeneracy of the magnetic spin ice^^. 
At least qualitatively, we have evidence that 

^inag arises 

from the monopole excitations and that its magnitude 
mainly depends on the monopole mobility which is re- 
lated to the degree of the ground-state degeneracy. How- 
ever, on a more quantitative level there are many open 
questions. Hence, a more fundamental theory is required. 



V. CONCLUSION 

In conclusion, we find that the low-temperature ther- 
mal conductivity of the spin-ice compound Dy2Ti207 has 
a pronounced magnetic contribution K,nag in zero mag- 
netic field. This Kmag has a strong magnetic-field de- 
pendence, which is highly anisotropic with respect to the 
magnetic-field direction. The magnitude of ^mag directly 
correlates with the degree of the ground-state degener- 
acy. It is maximum in zero field, where the ground-state 
degeneracy is maximum (6-fold). Applying finite mag- 
netic fields along different directions stabilize other (spin- 
ice) ground states with reduced degrees of degeneracy 
depending on the direction and magnitude of the mag- 
netic field. A 3-fold degenerate Kagome-ice phase occurs 
for B II [111] below 1 T, where Kmag is reduced compared 
to the zero-field value, but remains larger than Kmag in 
the 2-fold degenerate ground state which is realized for 
B\ \ [110]. Non-degenerate ground states are reached ei- 
ther for B II [001] or for B || [111] above 1 T and in both 
cases, Kmag vanishes completely. For B\ \ [110], we ob- 
serve an additional anisotropy of Kmag with respect to 
the direction of the heat current j . For a heat current di- 
rected along the a chains, whose spins are easily polarized 
because of their finite components along the field direc- 
tion, Kniag is suppressed more strongly as compared to a 
heat current j directed along the perpendicular oriented 
/3 chains, whose spins are _L B and are thus unaffected by 
the magnetic field. The anisotropy of Kmag with respect 
to the direction of j is, however, too weak to be explained 
within a single spin-flip formalism that only accounts for 
nearest-neighbor interaction. 

Independent of the direction of the magnetic field, 
there are pronounced hysteresis effects when the mag- 
netic field is reduced from larger values towards zero 
field in the low-temperature region. In all these cases, 
k{B — >■ 0) approaches a reduced value compared to 
obtained in the zero-field cooling process, and as a func- 
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tion of time, the reduced k{B — >• 0) value slowly relaxes 
towards Kq. This low-field hysteresis of k{B — s- 0) corre- 
lates with the occurrence of a small, but finite remnant 
magnetization observed for all studied field directions at 
low temperature. For however, we also ob- 

serve a hysteresis of Kmag that has no analogue in the 
magnetization data. Within the Kagome-ice phase, K,nag 
strongly depends on whether this phase is entered from 
zero or from high magnetic fields. Increasing the mag- 
netic field from the degenerate zero-field ground state re- 
sults in a Kmag of the Kagome-ice phase that is strongly 
reduced compared to the value of Kmag which is obtained 
when the Kagome-ice phase is entered from above, i.e. 
from the fully polarized state starting from higher fields. 
This larger Kmag is practically field-independent within 
the Kagome-ice phase and it neither depends on time, 
whereas the lower Kmag slowly relaxes towards the larger 
plateau value. 

Furthermore, we also identify a considerable field de- 
pendence of the phononic contribution Uph which emerges 



at higher field, i.e. above the (almost) saturation of 
the magnetization. Essentially, the same field-dependent 
Kph{B) is observed in the half-doped reference compound 
(Dyo.5Yo.5)2Ti207, where the spin-ice features are al- 
most completely suppressed and all the above-described 
low-field features of Kmag are essentially absent. The field 
dependence of Kph can be explained by lattice distortions 
due to torques on the non-collinear Dy momenta induced 
by the external magnetic field, and the resulting small 
tilts of the Dy momenta can also explain the continuous 
increase of the high-field magnetization. 
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